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Evolutionary constraints on the masses of the components 
of HDE 226868/Cyg X-1 binary system 
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ABSTRACT 

Calculations carried out to model the evolution of HDE 226868, under different as- 
sumptions about the stellar wind mass loss rate, provide robust limits on the present 
mass of the star. It has to be in the range 40 ± 5 Mq if the distance to the system 
is in the range 1.95 to 2.35 kpc and the effective temperature of HDE 226868 in the 
range 30000 to 31000 K. Extending the possible intervals of these parameters to 1.8 
to 2.35 kpc and 28000 to 32000 K, one gets for the mass of the star the range 40 
± 10 M0. Including into the analysis observational properties such as the profiles of 
the emission lines, rotational broadening of the absorption lines and the ellipsoidal 
light variations, one can estimate also the mass of the compact component. It has to 
be in the ranges 20 ± 5 Mq and 13.5 29 Mq for the cases described above. The 
same analysis (using the evolutionary models and the observational properties listed 
above) yields lower limit to the distance to the system of ~ 2.0 kpc, if the effective 
temperature of HDE 226868 is higher than 30000 K. This limit to the distance does 
not depend on any photometric or astrometric considerations. 

Key words: binaries: general - X-ray sources - evolution - stars: individual (Cyg 
X-1) - masses of the components. 
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1 INTRODUCTION 

Cyg X-1 was the first binary system in which the presence of 
a black hole weis suggested (Bolton 1972; Webster & Mur- 
din 1972). For about a decade, it was the only object of 
that type. Many, sometimes exotic, models and scenarios 
were devised to avoid the presence of a black hole and to re- 
place it with a neutron star. With the advent of subsequent 
black hole candidates (LMC X-3, A620— 00), the motivation 
for such efforts substantially diminished. At present, there 
are no doubts about the presence of a black hole in the 
system. However, in spite of three decades of the investi- 
gations, there is still substantial uncertainty concerning the 
masses of both components. The mass function is known 
rather precisely. Its most recent value was given by Gies et 
al. (2003): f{M^) = 0.251 ± 0.007 Mq. There are also no 
doubts that the optical component must be close to filling 
its Roche lobe. Gies & Bolton (1986a, 1986b), analyzing the 
emission lines of the stellar wind, the rotational broaden- 
ing of the absorption lines of the optical component and the 
photometric V-band light curve, found that the fill-out fac- 
tor must be greater than 0.9 and its best value is in the 
range 0.95 to l.The value of the inclination of the orbit is 
less certain. From the analysis mentioned above (rotational 
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broadening of the absorption lines of the optical compo- 
nent and the ellipsoidal light variations in the V-band light 
curve), the authors estimated the inclination of the orbit to 
be i = 33 ± 5°. On the other hand, from polarimetric mea- 
surements (in three colours) Dolan & Tapia (1989) found 
the inclination i — G2^^ °, which, while different, is not 
inconsistent with the value of Gies & Bolton. 

Paczyhski (1974) has shown that basing only on the 
lack of the X-ray eclipses and the fact that the optical com- 
ponent cannot be larger than its Roche lobe, one can obtain 
lower limits to the masses of both components as functions 
of the distance to the system. Paczyhski found that the mass 
of the compact component has to be larger than 3.2 Mq if 
the distance is greater than 1.3 kpc. He noticed also that, 
basing on stellar evolution consideration, the distance can- 
not be larger than about 2.1 kpc. The distance to Cyg X-1 
is not very well established, but, certainly, it is greater than 
1.3 kpc. Margon, Bowyer & Stone (1973) and Bregman et 
al. (1973) analyzed interstellar reddening in the direction of 
HDE 226868 and found that its distance is about 2.5 kpc. 
Wu et al. (1982), from a study of the UV colors of HDE 
226868, concluded that its distance must be greater than 
1.9 kpc. Bregman's value was confirmed by Ninkov, Walker 
& Young (1987), who got the distance of 2.5 ± 0.3 kpc from 
the equivalent width of the H^ line. HDE 226868 lies only 
1° from the center of NGC 6871 (the core of the Cyg 0B3 
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association) and seems to be its member. However, the es- 
timates of the distance to the association (or its core) gave 
mixed results. Crawford, Barnes & Warren (1974) used three 
methods (UBV photometry, the H-R diagram fitting and the 
equivalent width of the line) and obtained the distance 
to the association of about 2 kpc. Humphreys (1978) esti- 
mated the photometric distance to Cyg 0B3 to be about 
2.3 kpc. Janes & Adler (1982) quote in their catalogue the 
distance to NGC 6871 as only 1.8 kpc. The value of 1.8 kpc 
as the distance to Cyg OB3 was obtained again by Gar- 
many & Stencel (1992), who used the H-R diagram fitting 
method. The discrepancy between the estimates of the dis- 
tances to HDE 226868 and to Cyg OB3 cast the doubt on 
the membership of Cyg X-1 in NGC 6871. Fortunately, the 
extensive photometry and spectroscopy by Massey, Johnson 
& DeGioia-Eastwood (1995) seemed to clarify the situation. 
They estimated the distance to NGC 6871 as 2.14 ± 0.07 
kpc, which was in reasonable agreement with the indepen- 
dent (H-y width) estimate of the distance to HDE 226868 
(Ninkov ct al. 1987). Let us note that a very similar value 
of the distance (d « 2.15 kpc) was derived also by Gies & 
Bolton (1986a) as a byproduct of their careful analysis of 
a large set of observational data for HDE 226868/Cyg X- 
1 binary system. Their result (see the point Mx — 16 M0 
and Mopt = 33 M© in their fig. 10) was obtained quite inde- 
pendently of any photometric considerations. We should also 
note, however, that while Herrero et al. (1995) do not discuss 
explicitly the value of the distance in their atmosphere mod- 
elling paper, their final model implies the value ~ 1.8 kpc. 
Finally, the question of Cyg X-1 membership in NGC 6871 
was definitely solved in recent paper by Mirabel & Rodriguez 
(2003). Comparing the high precision VLBI astrometry for 
Cyg X-1 and Hipparcos astrometry for the members of Cyg 
0B3, the authors convincingly demonstrated that Cyg X- 
1 shares, quite precisely, common proper motion with Cyg 
0B3 association. The velocity of Cyg X-1 with respect to the 
Sun is 70 ± 3 km/s. The relative space velocity of Cyg X-1 
with respect to Cyg 0B3 is only 9 ± 2 km/s, which is typical 
of the random velocities of stars in expanding associations 
(Blaauw 1991). 

After the work by Massey et al. (1995), several more 
estimates of distances to both Cyg OB3 and HDE 226868 
were made and the results still exhibit a substantial scatter. 
Malysheva (1997), using different photometry than Massey 
et al. (1995), obtained for the distance to Cyg 0B3 again 
d ~ 1.8 kpc. Dambis, Mel'nik & Rastorguev (2001) calcu- 
lated the trigonometric distance to Cyg OB3 as a median 
value of the Hipparcos distances to the 18 individual mem- 
ber stars and found d = 2.3 (4-1.4,-0.6) kpc. In the case of 
Cyg 0B3 their trigonometric distance agreed very well with 
the photometric distance of Blaha & Humphreys (1989), al- 
though they found that, typically, Hipparcos trigonometric 
distances to OB associations are by ~ 12% smaller than pho- 
tometric distances of Blaha & Humphreys. A recent photo- 
metric distance estimate made by Mikolajewska (2002, pri- 
vate communication) gave the result d = 1.95 ±0.10 kpc. 
Lestrade et al. (1999) used the VLBI astrometry to esti- 
mate the distance to Cyg X-1 and obtained d = 1.0 2.3 
kpc. 

All this extensive effort leads to the conclusion that the 
distance to HDE 226868/Cyg X-1 binary system cannot be 
significantly different from 2 kpc. Therefore, I believe that 



the best choice will be to follow Massey et al. (1995) and use 
d = 2.15 ± 0.07 (1(7 error) or ± 0.2 (3a error) kpc as the 
distance to Cyg X-1. As we shall see at the end of the paper, 
this choice will be, quite independently, supported by evolu- 
tionary considerations. The above value is used throughout 
the rest of this paper (except when I assume the distance to 
be a free parameter) . In those parts of the discussion where 
the distance is assumed to be a free parameter, I shall con- 
sider the value 1.8 kpc as a reasonable lower limit to the 
distance of Cyg X-1. 



2 THE MODEL INDEPENDENT LOWER 
LIMITS TO THE COMPONENTS MASSES 

In this section, I shall repeat Paczynski's (1974) analysis, us- 
ing the present day data. This analysis is based only on solid 
observational facts and is, therefore, model independent. Let 
us remind that these solid facts include: (i) the value of the 
mass function, (ii) the lack of the X-ray eclipses, (iii) the ob- 
vious fact that the optical component cannot be larger than 
its Roche lobe, and (iiii) the spectral type and the photom- 
etry of the optical component (used to derive its effective 
temperature and the bolometric correction). 

The most recent determination of the mass function is 
f{M^) = 0.251 ± 0.007 Mq (Gies et al. 2003), which differs 
only slightly from the earlier values f{Mx) = 0.244 ± 0.005 
Mq (Brocksopp et al. 1999) and /(M^) = 0.252 M© (Gies 
& Bolton 1982). The corresponding projected radius of the 
optical component orbit is aisini = 8.36 ± 0.08 R©. The 
spectral type of HDE 226868 was determined by Walborn 
(1973) as 09.7 lab. This classification was confirmed by Gies 
& Bolton (1986a), by Ninkov et al. (1987) and by Herrero 
et al. (1995) and was used in most of the literature. On the 
other hand, Massey et al. (1995) obtained, from their new 
massive spectroscopy of OB associations, the somewhat ear- 
lier spectral type - ON9 Ifa-|-. According to the calibration 
of effective temperatures for late O-type supcrgiants derived 
by Vacca, Garmany & ShuU (1996, hereafter VGS), one has 
Te = 32740 K for 09 la star and extrapolation for 09.7 
la gives Tc = 30690 K. On the other hand, Herrero et al. 
(1995) used atmosphere models to reproduce the observed 
spectrum of HDE 226868 and obtained Te = 32000 K as the 
best fit for the effective temperature of this star. 

Both VGS and Herrero et al. calibrations were based 
on pure H-He models of the atmospheres that neglected the 
effects of metallic lines blanketing. It is well known by now, 
that taking into account the line blanketing leads to the 
lower effective temperature for a given spectral type (Mar- 
tins et al. 2002, Repolust et al. 2004, Markova et al. 2004). 
The most recent calibration of effective temperatures for O- 
type supergiants (accounting for the line blanketing effects) 
was derived by Markova et al. (2004) . Their scale gives gen- 
erally lower effective temperatures than VGS scale (by up 
to 10 000 K for the earliest spectral types). However, for 
the spectral type 09.7 I, both temperature scales converge 
to the value Te = 30700 K. It seems that, at the present 
state of knowledge, the most reasonable approach would be 
to assume that To = 30700 K is the best estimate of the 
effective temperature of HDE 226868. To take into account 
the uncertainty of this estimate, in further discussion I shall 
consider a broad interval of Te = 28000 ^ 32000 K, as the 
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range of the possible values of the effective temperature of 
HDE 226868. 

As the effective temperature-bolometric correction re- 
lation is not affected by the line blanketing (Martins et 
al. 2002), it is possible to take this relation from VGS 
tables. For the effective temperature 30700 K, the appro- 
priate value of the bolometric correction is B.C. = —3.06 
. The corresponding unreddened colour index should be 
{B - V)o = -0.26 (Schmidt-Kaler 1982). The most recent 
photometry of HDE 226868 by Masscy ct al. (1995) gives 
V = 8.81 and B - V = 0.83. Therefore the reddening is 
Eb-v = 1.09 and the unreddened V magnitude is Vo = 5.43 
(following Masscy ct al., I adopt _Rv ~ 3.1). With the 
effective temperature and the bolometric correction given 
above, wc can then express the radius and luminosity of 
HDE 226868 as functions of the distance: 

R^pt/Re = 10.59 (d/lkpc), (1) 

Lopt/LQ = 8.95 X 10* (d/lkpc)^. (2) 

Following the procedure of Paczynski (1974), we can 
now calculate the lower limits to the masses of both compo- 
nents as functions of the distance to the system (under the 
assumptions listed at the beginning of this section). The re- 
sults of these calculations are given in the first part of Table 
1. 

Let us notice that the changes, comparing with 
Paczyhski's Table 1, arc mostly due to substantial change 
in the effective temperature (Paczynski took Te = 25000 
K) and the corresponding change in the bolometric correc- 
tion for HDE 226868. Let us also notice that the effects of 
the changes in the effective temperature and the bolometric 
correction, to a large degree, cancel each other when we con- 
sider the radius of HDE 226868 and the lower limits to the 
masses of both components. They introduce, however, the 
dramatic (by a factor of ^ 2.5) increase of the luminosity of 
HDE 226868 (for a given distance). As a resuh, HDE 226868 
appears to be a very bright star: its optical luminosity is ~ 
1 2 X 10^® erg/sec (by two orders of magnitude higher 
than the typical X-ray luminosity of Cyg X-1). 

Undoubtedly, we know now the calibration of effective 
temperatures for O-type supergiants much better than 30 
years ago. However, still, there remains some uncertainty, 
clearly demonstrated by the recent major revision intro- 
duced by accounting for the line blanketing effects (as dis- 
cussed earlier). To investigate the effect of this uncertainty, 
I assumed that the effective temperature of HDE 226868 
lies in the range 28000 32000 K and calculated the cor- 
responding lower limits to the masses of both components 
(using the modified versions of eqs. (1) — (2)). The results of 
these calculations are given in the second and the third part 
of Table 1. All sets of results are illustrated in Fig. 1. 

Comparing the both sections of the table (or looking 
at the Fig. 1), we may notice that the uncertainty of the 
effective temperature does not affect significantly the lower 
limits to the masses. 

Summarizing the considerations of this section, we may 
state that the masses of the compact and the optical com- 
ponents must be greater than ^ 8 Mq and ~ 29 Mg, re- 
spectively (if the distance to the system is in the range 2.15 
± 0.07 kpc). Let us remind that these values are model in- 
dependent and therefore very hard to contest. If one wants 



Table 1. Parameters of the components of HDE 
226868/Cyg X-1 binary system as functions of the as- 
sumed distance d for three values of the effective temper- 
ature of the optical component. 
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to be more conservative and uses the 3a error as the un- 
certainty of the distance estimate, then these lower limits 
become ~ 7 Mq and ~ 24 M©, respectively. Even adopting 
the smallest value ever claimed for the distance (d = 1.8 
kpc), one still obtains the limits ~ 6 Mq and ~ 19 M©, 
respectively. One should add that our method gives only 
approximate estimates, since (as noted by Gies 2004, pri- 
vate communication) it assumes isotropic optical emission, 
which is not the case for tidally distorted star. The resulting 
inaccuracies should not be, however, significant. 

The lower limits, presented in Table 1 and Fig. 1, be- 
come the actual values if the fill-out factor is equal 1.0 and 
the inclination of the orbit corresponds to the grazing eclipse 
orientation. The fill-out factor cannot be much different from 
1.0 (Gies & Bolton 1986a,b). However, the inclination is the 
source of larger uncertainty. If the inclination is close to 62° 
(as suggested by the polarimetry, Dolan & Tapia 1989), then 
we are close to the grazing eclipse situation and the masses 
are close to the lower limits obtained in this section, i.e. 8-^9 
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Figure 1. The lower limits to the masses of the optical (A) and 
compact (B) components as functions of the assumed distance, 
d. Thick solid lines correspond to the most likely value of the ef- 
fective temperature of HDE 226868 (Tc = 30 700 K). The broken 
lines and the dash-dotted lines correspond to the effective tem- 
peratures 28 000 K and 32 000 K, respectively. Thin solid vertical 
line indicates the most likely value of the distance {d = 2.15 kpc). 
The broken vertical lines indicate ±lcr errors in the distance es- 
timate. The ±3(T error range corresponds to the distance interval 
1.95 2.35 kpc. 

Mq and 30 36 M0 , respectively. If (as it seems more likely 
— see the discussion in Section 5), the inclination is rather 
close to 33° (advocated by Gies & Bolton 1986a and still not 
excluded by polarimetry) , then the masses are substantially 
higher and probably close to the values suggested by Gies 
& Bolton (16 ± 5 Mq and 33 ± 9 Mq, respectively). I shaU 
return to this point in further discussion. 



3 THE EVOLUTIONARY STATUS OF HDE 
226868 

Quite independent constraints on the masses may be ob- 
tained from analysis of the evolutionary status of the op- 
tical component. First, let us note that we know relatively 
well its luminosity. At the distance of 2.15 ± 0.07 kpc, the 
luminosity must be Mboi « —9.29 ± 0.07 (equivalent to 
log(L/Lo) = 5.62 ± 0.03 in Table 1). Including the uncer- 



tainty of the effective temperature and assuming 3o" error in 
the distance estimate, the range for the luminosity becomes 
Mboi « -8.8^-9.6 (equivalent to log(L/LQ) = 5.40 5.75 
in Table 1). The real range is somewhat wider because one 
should include also the errors in Eb-v (about ± 0.05) and 
Ry (about ± 0.1) estimates. The final range becomes there- 
fore Mboi w —8.5 -. — 9.9. One may note that this range 
corresponds to My w —5.7 -= — 6.7, which agrees quite well 
with Mv « -6.5±0.2 obtained by Ninkov et al. (1987) from 
the equivalent width of the line. 

Second, let us consider the evolutionary phase of HDE 
226868. Even without the specific evolutionary calculations, 
it is relatively straightforward to argue that it must be a core 
hydrogen burning configuration. It comes from the fact that 
HDE 226868 must be in a relatively stable phase of its evo- 
lution. We observe no measurable variations of the orbital 
period of the system. The 3cr upper limit for the relative 
change of the orbital period |d In P/d f | is ~ 2 3 x 10~^ 
yr~^ (Gies & Bolton 1982). This implies that the present 
evolutionary timescale of HDE 226868 must be substantially 
longer than ~ 3 5 x 10"* yr. This excludes the possibility 
of the post-main sequence evolution, as the expected evo- 
lutionary time scale is, in that case, of the order of 10^ yr. 
Ziolkowski (1977), investigating the post- main sequence evo- 
lution of a 25 Mq star (of similar radius but smaller lumi- 
nosity than HDE 226868) approaching its Roche lobe, found 
that its radius was increasing from the fill-out factor equal 
0.95 to 1.00 in only ~ 600 yr. After filling the Roche lobe, it 
took only ~ 200 yr for the mass transfer rate to exceed the 
Eddington limit by 3 orders of magnitude. For HDE 226868, 
which is more luminous and, probably, more massive (35-^45 
M0, as we shall see later), the relevant time scales should be 
even shorter. If HDE 226868 were expanding that fast (5% 
radius increase in less than 600 yr), we should observe the 
noticeable rise of its stellar wind strength and of the X-ray 
luminosity of the system (both parameters are sensitive to 
the fill-out factor if it is close to 1). Moreover, the chance 
of observing the system during such a short time window 
(just few hundred years, before it gets extinguished as an X- 
ray source due to hyper-Eddington accretion) is very small. 
Therefore, we have to assume that HDE 226868 (similarly 
as the optical components of other massive X-ray binaries - 
see Ziolkowski 1977) must still burn hydrogeu iu its core. 

Exactly the same conclusion may be obtained directly 
as a unique outcome of the numerical modelling of the evo- 
lution of HDE 226868. This topic will be discussed in the 
next section. 



4 THE EVOLUTIONARY CALCULATIONS 
FOR HDE 226868 

4.1 The general description 

I computed evolutionary tracks for core hydrogen burn- 
ing phase of stars with the initial masses in the range 
40 -i- 80 M0. The Warsaw evolutionary code developed by 
Bohdan Paczyhski and Maciek Kozlowski and kept updated 
by Ryszard Sienkiewicz was used. The initial chemical com- 
position X=0.7 and .^=0.03 is adopted. The opacity tables 
incorporating OPAL opacities (Iglesias & Rogers 1996) as 
well as molecular and grain opacities (Alexander & Fergu- 
son 1994) were used. The nuclear reaction rates are those 
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of Bahcall & Pinsonncault (1995). The equation of state 
used was Livermore Laboratory OPAL (Rogers, Svenson & 
Iglesias 1996). I neglected the semiconvective mixing, as it 
is not important during the evolutionary phase considered 
(most of the models of interest had central hydrogen content 
Xc > 0.2). Similarly any overshooting at the border of the 
convective core was neglected (it is even less important). 

The calculations were carried out under the assumption 
that the evolution starts from the homogeneous configura- 
tions. It means that the consequences of the fact that some 
of the matter of the star, possibly dumped from the progen- 
itor of the present black hole, could have somewhat mod- 
ified chemical composition, were neglected. It means also, 
that the consequences of the fact that some nuclear evo- 
lution (hydrogen burning) could, possibly, take place while 
the mass of the star was smaller (prior to the mass transfer) 
were neglected as well. It seems that both simplifications 
do not alter significantly the outcome of the evolutionary 
calculations. I shall return to this point in later discussion. 



4.2 The stellar wind mass loss 

The most uncertain clement of the calculations of the early 
evolution of massive stars is the mass loss due to stellar 
wind. The uncertainty of the estimate of its rate is the sin- 
gle most important factor influencing the outcome of the 
calculations. The observations seem to indicate that there is 
a substantial scatter of the mass loss rates (typically, by a 
factor of two, but sometimes this factor can reach up to five) 
among the stars of similar luminosities and effective temper- 
atures (see e.g. de Jager, Nieuwenhuijzen & Van Der Hucht 
1988). I decided to use the prescription given by Nieuwen- 
huijzen & de Jager (1990) in the form of the fornnila derived 
by Hurley, Pols & Tout (1999, hereafter HPT). Bearing in 
mind that the formula gives the mass loss rate estimate with 
the accuracy that is probably not better than within a fac- 
tor of two, I introduced the multiplying factor /sw applied 
to HPT formula. For each initial mass of the star three evo- 
lutionary sequences were calculated with the value of the 
parameter /sw equal, in sequence, 0.5, 1 and 2. In this way, 
the uncertainty of the theory of evolution could be, hope- 
fully, taken into account. It might be interesting to note, at 
this point, the astonishingly good agreement between the 
HPT formula and the observational determination of the 
mass loss rate for HDE 226868: for the most likely values of 
the parameters of the star, the HPT value agrees with the 
observed one {M = -2.6 x 10"^' Mo/yr, Gies et al. 2003, 
see below) to an accuracy of about 5% (this is, definitely, 
better than the precision of both the HPT formula and of 
the observational estimate). 

I should remind at this point that the supergiant com- 
ponents in some high mass X-ray binaries are significantly 
undermassive for their luminosities (Ziolkowski 1977). In 
some systems, like Cen X-3, this undermassivness is very 
serious and requires much stronger mass loss than the nor- 
mal stellar wind (Ziolkowski 1978). However, most likely, it 
is not the case for HDE 226868. Its parameters may be fully 
explained by the evolution with the normal stellar wind (as 
will be demonstrated in the further discussion). 



4.3 The evolutionary tracks 

Some of the obtained evolutionary tracks in the H-R dia- 
gram are shown in Fig. 2. The careful reader might notice 
that the luminosities of my models are by about 0.1 dex (~ 
25%) smaller than the evolutionary tracks of Schalller al. 
(1992). This difference should be attributed mainly to the 
fact that I used later edition of the opacity tables contain- 
ing higher values of the opacities (the opacity tables always 
evolve in the direction of the growing opacity - never the 
other way). The correctness of the new opacities was con- 
firmed by stellar pulsations calculations - see e.g. Pamyat- 
nykh (1999). Part of the luminosity difference is due to the 
fact, that I used higher metallicity {Z = 0.03 instead of Z= 
0.02). And, finally, part of the luminosity difference (all fac- 
tors work in the same direction) is due to the fact, that my 
stellar winds are stronger (again due to higher metallicity). 

As may be seen from a quick look at Fig. 2, the initial 
(zero age main sequence) mass of HDE 226868 had to be 
in the range 35 to 55 Mq. The masses of the models cor- 
responding to the present day state of HDE 226868 (inside 
the 3(T parallelogram) are in the range 32.4 to 50.5 M© . The 
central hydrogen content in these models is between 0.126 
and 0.265 and their evolutionary age (since the beginning of 
the central hydrogen burning) is between 2.7 and 4.3 xlO® 

yr- 

To obtain evolutionary models which satisfactorily re- 
produce the present day state of HDE 226868, one has to 
match not only the luminosity and the effective tempera- 
ture but also the rate of the stellar wind mass loss. The 
first observational determination of this parameter for HDE 
226868 was done by Hutchings (1976) who, analyzing the 
visual spectrographic data, got the value M — —2.5 x 10~® 
Mo/yr. Persi et al. (1980) estimated the rate of the mass 
outflow from the infrared emission of the expanding circum- 
stellar envelope and got the value M = —3.5 x 10^'' Mq/jt. 
More recent estimate of Herrero et al. (1995) was based on 
the fits of the Ha profile and led to the conclusion that the 
mass loss rate Ues between 2 and 6 xlO~® Mo/yr. They 
noted that the inaccuracy of the fits is, probably, due to 
the fact that the stellar wind from HDE 226868 is focused 
towards the black hole (for the discussion of the "focused" 
wind model see Gies & Bolton 1986 and Miller et al. 2002). 
In their test modelling of the atmosphere of HDE 226868, 
Herrero et al. were using the value M = —4 x 10~® Mo/yr 
and for their final model they chose M — —3 x 10~^ M© /yr. 
The most recent estimate by Gies et al. (2003) gave the value 
M = —2.6 X 10~® Mo/yr for the low/hard state (which is a 
typical state of Cyg X-1) and M = -2.0 x 10"*^ Mo/yr for 
the high/soft state (which is less frequent in this source). 
Taking all this into account, I assumed that the observed 
rate of mass outflow from HDE 226868 is M = -2.6 x 10"" 
Mo/yr with an error by a factor of about 2, i.e. I assumed 
that the observed rate of mass outflow is between 1.3 and 
5.2 x lO^*" Mo /yr. Let me remind that our theoretical evo- 
lutionary models were using mass loss rates in the range 0.5 
to 2 times the rate given by HPT. Altogether, it means, that 
the model for which the mass loss rate calculated with the 
HPT formula would be by a factor of up to four smaller or by 
a factor of up to four larger than the nominal observational 
value (M = -2.6 x 10~^ Mo/yr) is still considered to be 
an acceptable match. This is, probably, more than sufficient 
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Figure 2. The evolutionary tracks in the H-R diagram. The traeks are labeled with the initial mass of the star (in solar units). The solid 
lines describe the tracks computed with the stellar wind mass loss rates according to HPT (Hurley et al., 2000) formula. The broken lines 
and the dash-dotted lines describe the tracks computed with the mass loss rates smaller by a factor of two and larger by a factor of two, 
respectively. The slanted dotted lines correspond to the position of HDE 226868 for different assumed values of its distance (the assumed 
value of the distance in kpc is given at the right end of each line). The vertical dotted lines correspond to the effective temperatures 
of HDE 226868 equal (from left to right) to 32, 31, 30 and 28 xlO^K. The most likely position of HDE 226868 lies within the large 
parallelogram (±3(T error in distance). 



allowance for the uncertainty of the mass loss rates. The pa- 
rameters of some of these models are given in Table 2. The 

positions of the acceptable models in the mass— luminosity 
diagram (for the most likely value of the effective tempera- 
ture Te = 30700 K) are shown in the Fig. 3. Since the lu- 
minosity of a model of HDE 226868 can be (for an assumed 
effective temperature) directly translated into the distance 
to the system (see Table 1), the Fig. 3 can be also considered 
as the mciss— distance diagram. The corresponding distance 
scale is shown on this picture. The line corresponding to 
the lower mass limit— distance relation obtained in section 
2 (see Table 1 and Fig. 1) is also shown. One may notice 
that the deduced mass of HDE 226868 depends mainly on 
the assumed distance to the binary system, and only very 
weakly on the assumptions about the mass loss rates. For 
the assumed effective temperature (Tc = 30700 K) this mass 
has to be in the range 35 -^ 45 M© (if the distance to the sys- 
tem is in the range 1.8 to 2.35 kpc). One may also conclude 
(with the help of the lower mass limits derived in Section 2) 



that the distance to the system cannot be larger than ~ 2.5 
kpc (still for the assumed effective temperature). 

Similar diagrams as Fig. 3 may be constructed for other 
possible values of the effective temperature of HDE 226868. 
Qualitatively, they look similar to the Fig. 3. For the effec- 
tive temperature Tc = 28000 K the mass of the star has to 
be in the range 29 ^ 35 Mq and the upper limit for the 
distance to the system is ~ 2.1 kpc. For the effective tem- 
perature Te = 32000 K the mass of the star has to be in the 
range 37 50 M0 and the upper limit for the distance to 
the system is ~ 2.7 kpc. 

5 DISCUSSION 

5.1 The mass of HDE 226868 

The obvious and fairly strong conclusion derived from the 
evolutionary calculations is that HDE 2226868 had to be 
quite massive in the past and is still very massive at present. 
To put it very briefly, it has to be very massive because it is 
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Figure 3. The positions of the evolutionary models of HDE 226868 (with an assumed effective temperature Te = 30700 K) in the 
mass— luminosity and mass— distance diagram ("mass" means here the present mass of the star). The scale of the distances (in kpc) 
is shown with the thin broken vertical lines (the most likely value of the distance is 2.15 ± 0.2 kpc, 3o" error). The triangles, circles 
and squares correspond to the models with the stellar wind mass loss rates calculated with the HPT formula multiplied by the factors 
0.5, 1 and 2, respectively. Only the models with the acceptable rates of mass loss (in the range 0.25 4 times the nominal observed 
value, M = —2.6 x 10~® MQ/yr) are shown. The solid line (without circles) describes the relation between the distance and the lower 
limit for the mass of HDE 226868 obtained in Section 2. The location of crossing of this relation by the sequences of the evolutionary 
models indicates that the distance to the binary system cannot be larger than ~ 2.5 kpc (for the assumed effective temperature). The 
oval indicates the range of the models of the optical star for which viable models of the binary system could be constructed (see the 
discussion in the Section 5.2). 



very bright (while still burning hydrogen in the core). For the 
most likely range of parameters (30000 to 31000 K for the 
eflFective temperature and 2.08 to 2.22 kpc for the distance) 
the present mass of HDE 2268668 has to be in the range 
37 to 44 Mq. Increasing the interval of the possible values 
of the distance to 1.95 to 2.35 kpc increases this range only 
slightly (to 33.5 -^47Mq). Extending the interval of the pos- 
sible values of the effective temperatures to 28000 to 32000 
K leads to the range of possible masses 31 -tSOMq. Finally, 
extending the range of the possible distances down to 1.8 
kpc (the smallest value quoted in the literature) results in 
the final range of the possible masses of HDE 226868: 29 
^50Mq . The value of 29 M© seems to be a firm lower limit 
for the present mass of HDE 226868. 

One may ask whether it is possible to construct substan- 
tially less massive evolutionary models by taking more lib- 
eral estimate of the observational uncertainties or account- 



ing for some simplifications of our models. The answer is neg- 
ative. Decreasing the effective temperature of HDE 226868 
would result in lower luminosities and so would lead to less 
meissive configurations. However, the value Te = 28000 K, 
taken as the lower limit to the effective temperature in the 
above consideration is probably already too low (the true 
value is probably in the range 30000 to 31000 K). Decreas- 
ing the minimum distance from 1.95 to 1.8 kpc decreased the 
minimum possible masses of the models by only ~ 2M0, as 
was mentioned above. Let us now consider the possible con- 
sequences of some simplifications of our evolutionary cal- 
culations. I neglected the possible dumping of matter on 
HDE 226868 from the progenitor of Cyg X-1 (and so "re- 
juvenation" of the star) during the first mass loss/exchange 
phase. Most likely, there was no significant accretion on HDE 
226868 because any substantial mass transfer in such a close 
binary and in so early stage of evolution (when both stars 
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had no distinct cores) would lead to the rapid built-up of a 
common envelope and the coalescence of both components. 
So, most likely, the first mass transfer was mainly the mass 
loss from the system. If there were no serious accretion, then 
the nuclear evolution prior to the mass transfer is prop- 
erly accounted for in our calculations. If (which is rather 
unlikely) there was a substantial mass dumping on HDE 
226868, then the nuclear evolution before the mass trans- 
fer is negligible because the substantially less massive star 
was evolving much more slowly (large mass gain essentially 
resets the evolutionary clock to ZAMS). 

One may wonder whether smaller metallicity {Z = 0.02 
instead of Z= 0.03) would not be a better choice for HDE 
226868. To check the possible effects of such alternative 
choice, I calculated several evolutionary tracks for stars with 
Z = 0.02. As might be expected, the changes in the results 
were small. For the models reproducing the "best" parame- 
ters of HDE 226868, the masses decreased from about 41.7 
Mo to about 39.9 Mq (only about 4%). 

There remains a rather remote possibility that the 
chemical composition of HDE 226868 is not normal, i.e. that 
it contains less hydrogen than a normal Population I star. 
E.g. Horrero et al. (1995) concluded, from fitting models 
of the atmosphere, that helium might be overabundant by 
a factor of about two. The evolutionary calculations pre- 
sented here indicate that the stellar wind mass loss from 
HDE 226868 was not severe enough for hydrogen depleted 
layers to show up on the surface. However, if a large amount 
of matter with significantly decreased hydrogen content was 
dumped on HDE 226868 from the progenitor of Cyg X-1, 
then the star might contain less hydrogen than the assumed 
value of 0.7. The lack of the CNO anomalies (Dearborn 1977) 
testifies against this being the case. However, to check the 
consequences of such (rather unlikely) situation, I calculated 
several evolutionary sequences for the stars with the initial 
chemical composition X—0.5, Z—0.03 (this corresponds to 
an overabundance of helium by a factor of two). As could 
be easily expected, the models with similar luminosities were 
now less massive, but the change was not very substantial. 
For the "best fit" models, the masses decreased from about 
41.7 Mc5 to about 35.3 M., (only about 15%). However, I 
would like to stress again, that there are no good reasons to 
expect such abnormal hydrogen content in HDE 226868. 

The only way to produce an evolutionary model with 
the mass < 20 Mq is to decrease drastically the distance to 
the system. I constructed one evolutionary track that pro- 
duced a 19.5 M0 configuration at the effective temperature 
30700 K with the mass loss rate M « -2.5 x 10"® (the ini- 
tial mass of the star was 32.5 M©). However, the price was 
very high: the distance had to be decreased to ~ 0.9 kpc 
and the rate of the mass loss during the evolution had to be 
multiplied by a factor of ~ 20 with respect to the HPT for- 
mula. Any of these conditions, even taken separately, seems 
extremely unlikely. 

There is a clear confiict between the above considera- 
tions and the mass estimate given by Herrero et al. (1995). 
They used atmosphere models to reproduce the observed 
spectrum of HDE 226868. Then, they used their fit to de- 
termine the effective temperature, the surface gravity, the 
radius and the helium abundance of the star. Their mass 
estimate (M ~ 18 ± 4 M©) is a consequence of the surface 
gravity and the radius determination. 



I see no solution of this conflict. I may only repeat that 
the evolutionary calculations for such an early evolution- 
ary phase are very robust and rather difficult to dispute 
(it is, essentially, the mass-luminosity relation for the main 
sequence stars). I may also add that the evolutionary re- 
sults presented above are roughly consistent with the model 
of Gies & Bolton (1986a), based on an extensive analysis of 
the large and diversified collection of the observational data. 

5.2 The Mass of Cyg X-1 

In the previous section we concluded that the present mass 
of HDE 226868 has to be in the range of 29 ^ 50 Mq (assum- 
ing the widest possible ranges of the distances and the ef- 
fective tempeatures). The mass of its companion, black hole 
Cyg X-1, does not come out directly from the evolutionary 
calculations. However, it can be calculated (with the help 
of some observational constraints), once we select a chosen 
evolutionary model of HDE 226868. Once the selection is 
made, we know the mass Mopt and the radius Ropt of the 
optical component. Wc can also calculate the distance, with 
the help of cq. (1) or (2) (or the corresponding expressions 
for other effective temperatures). Subsequently, we can use 
two equations to solve for the inclination of the orbit i and 
the mass ratio q = Mopt/Mx. One of these equations makes 
use of the mass function, 

/(Mx) = Moptsin"i/[a(l + if]- (3) 

The other relates the radius of the star to the size of the 
orbit, 

Ropt = RrL X /rl = 

= /RL(0.38-|-0.21og g)A = 

= /RL(0.38-h0.21og 9)oi(l-h9), (4) 

where _Rrl is the radius of the Roche lobe around HDE 
226868, /rl is the fill-out factor (/rl = Ropt/RRh), A is 
the orbital separation of the binary components and ai is 
the radius of the orbit of HDE 226868 (sec section 2). 

Inserting the observational data, eqs. (3) — (4) can be 
written as: 

Moptsin'i/[9(l + (if] = 0-251, (5) 

Ropt = /RL(0.38-H0.21og q){l + q) x 8.36/sini. (6) 

Once a given evolutionary model of HDE 226868 is se- 
lected from the grid of the acceptable models and a value 
of the parameter /rl is assumed, the eqs. (5) — (6) can be 
solved for i and q. Knowing q we can immediately calcu- 
late also the mass of compact component A/x- In principle, 
this procedure can be applied to any combination of the 
evolutionary model and of the value of /rl- In fact, how- 
ever, not every evolutionary model of HDE 226868 (accept- 
able if we consider the optical component alone) permits 
the construction of a consistent model of the binary system. 
This is because of the observational constraints on the value 
of the inclination i and, especially, because of the strong 
observational constraints on the value of the fill-out factor 
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/rl- As demonstrated by Gies & Bolton ( 1986a, b), in order 
to explain quantitatively the He emission lines produced in 
the stellar wind from HDE 226868, the fill-out factor /rl 
has to be larger than 0.9 and, most likely, not smaller than 
0.95 (perhaps the best value would bo around 0.98). On the 
other hand, Gies & Bolton demonstrated that observed ro- 
tational broadening of the photospheric absorption lines of 
HDE 226868 and the observed amplitude of the ellipsoidal 
light variations impose substantial constraints on the val- 
ues of the mass ratio, the fill-out factor and the inclination. 
The rotational broadening depends mainly on the mass ratio 
and the measured broadening indicates the mass ratio in the 
range ^ 2 ^ 2.5. The amplitude of the ellipsoidal light vari- 
ations is determined mainly by the fill-out factor and the 
inclination. For the assumed values of /rl equal 0.9, 0.95 
and 1, the resulting inclination is ~ 38°, ~ 33° and ~ 28°, 
respectively. 

Constructing models of the binary system to be consis- 
tent with observations, I assumed that for any adopted value 
of /rl, the calculated inclination should be within ±5° from 
the corresponding values quoted above. I started with dif- 
ferent models of the optical component, acceptable from the 
point of view of the stellar evolution, as described in sec- 
tion 4.3. Then, I assumed the value of /rl equal 0.95 and 
solved eqs. (5) — (6) to find q, i and Mx. Subsequently, I tried 
higher values of /rl- The higher values of /rl produced the 
solutions with lower (in many cases unacceptably low) val- 
ues of the inclination. The dependence of q, i and Afx on 
the assumed value of /rl may be seen from many sequences 
of binary models, presented in Tab. 2. For all sequences in- 
cluded in the table, I present only two limiting models for 
the lowest and the highest value of /rl for which an ac- 
ceptable model could still bo obtained. If there is only one 
entry for a given sequence {as is the case for Mopt = 39.67 
M© and Mopt = 47.37 M©), it means that only one value of 
/rl produced an acceptable solution (the value of /rl was 
varied with the increment of 0.01). 

I classified, as acceptable, the models satisfying the fol- 
lowing criteria: 

(1) d= 1.8 -=-2.35 kpc 

(2) To = 28000 ^ 32000 K 

(3) -M = 1.3 ^ 5.2 X 10"® Mo/yr 

(4) /rl > 0.95 

(5) i = 28° -h (1 - /rl) X 100° ± 5° (this corre- 
sponds to i « 28 ^ 38° for /rl = 0.95 and i 23 ^ 33° for 
/rl = 1.00, as advocated by Gies & Bolton 1986a). 

Parameters of selected acceptable models are given in 
the second part of Table 2. For each assumed value of the 
effective temperature, I present two sequences of models cor- 
responding to the lowest and the highest current mass of the 
optical component. The exception to this rule is the sequence 
of models for = 30700 K and Mopt = 41.67 M© (d = 2.15 
kpc). The initial evolutionary parameters (Mo and /sw) of 
the optical component model for this sequence were adjusted 
so as to obtain the perfect fit with the "best" values of the 
observational parameters of HDE 226868 (effective temper- 
ature, luminosity and the rate of stellar wind mass loss). 
All models from this sequence (for /rl = 0.96 1.00) are 
successful from the point of view of our criteria. 

There are several conclusions that can be drawn from 
the collection of the obtained models of the binary system 
(only some of these models are shown in Table 2). The first 



concerns the mass of the compact component. Assuming the 
widest possible intervals of the distances and the effective 
temperatures: 1.8 to 2.35 kpc and 28000 to 32000 K, the 
mass of the compact component must be in the range 13.5 
-^ 28.5 M©. For the most likely intervals of these param- 
eters: 1.95 to 2.35 kpc and 30000 to 31000 K, this range 
narrows to 15.5 25 M©. The second conclusion is related 
to the distance to the binary system. It appears that for 
a given (assumed) effective temperature, consistent mod- 
els are possible only for some, relatively narrow, interval of 
distances (the distances out of this interval would require 
unacceptably low or unacceptably high values of the incli- 
nation). For the effective temperatures equal to 28, 30, 30.7, 
31 and 32 xlO'^K, the corresponding intervals arc 1.80^1.88, 
1.98^2.15, 2.04^2.22, 2.07-=-2.27 and 2.21^2.35 kpc, respec- 
tively. It is certainly encouraging that the most likely range 
of the effective temperatures (30000 to 31000 K) requires the 
distance interval (1.98-^2.27 kpc) that almost exactly coin- 
cides with the independent estimate of the most likely dis- 
tance range (1.95-^2.35 kpc). The obtained relation between 
the distance and the effective temperature of HDE 226868 
means that, if in the future the distance will be known more 
precisely (e.g. from the future astromctric space missions), 
then it will be possible to set constraints on the effective 
temperature. For example, if it is found that the distance 
to the system is greater than 2.1 kpc, it would mean that 
the effective temperature of HDE 226868 has to be higher 
than 30000 K. Also the opposite is true. If it is found (e.g. 
from the better models of the atmospheres) that the effec- 
tive temperature of HDE 22868 is greater than 30000 K, it 
would mean that the distance to the system must be larger 
than 2.0 kpc (and this conclusion would not depend on any 
photometric or astromctric considerations). 

Finally, the third conclusion confirms the earlier results 
(based on the evolutionary calculations alone) limiting the 
present mass of HDE 226868 to the range 29 50 M© and 
the initial mass to the range 33.5 -^ 55 M©. 

Let us note that for the black hole mass in the range 
15 25 M©, the state transitions occurs in Cyg X-1 at the 
luminosity level equal 0.025 to 0.015 of the Eddington lumi- 
nosity (assuming d = 2.15 kpc and using the fiux values given 
by Zdziarski et al. 2002). These values lie in the range 0.007 
to 0.03 found for other X-ray binaries (Maccarone 2003). 

To summarize this section, the evolutionary consider- 
ations provide much stronger lower limits to the masses of 
both components than the model independent analysis, pre- 
sented in section 2. In addition, they yield some independent 
constraints on the distance to the system. 



6 THE POSSIBLE EVOLUTIONARY 
SCENARIO 

Let us speculate a little bit about the possible evolutionary 
past of our binary system. It had to start as a very massive 
system (the initial primary had to complete its evolution by 
the time the secondary (HDE 226868) reached its present 
evolutionary state i.e. in less than 3 -t- 4 x 10® yr) . The 
initial masses of the components were probably ^ 80 -i- 100 
and ~ 40-^50 M©. The massive primary (the progenitor of a 
black hole) was shedding the mass mainly in the form of the 
stellar wind (there was probably little, if any, mass transfer 
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to the companion). The mass of the primary decreased, at 
the end of its evolution (before the collapse), to ~ 15 -7- 25 
Mq. The low spatial velocity of Cyg X-1 with respect to its 
parental association indicates, as argued by Mirabel & Ro- 
driguez (2003), that the final collapse to a black hole pro- 
ceeded with very little (< 1 M0) mass ejection. It could be, 
even, a prompt collapse with no accompanying supernova 
explosion at all ("formation of a black hole in the dark"). 
The secondary (the progenitor of HDE 226868) was also los- 
ing mass in the form of the stellar wind, although at a lower 
rate. This mass loss took place both before and after the 
collapse of its companion. This process decreased the mass 
of the secondary to the present value of ~ 35 ^ 45 Mq . In 
this way, the present day binary consisting of a massive su- 
pergiant of ~ 40 ± 5 Mq (with a central hydrogen content 
of ~ 0.2 0.3) and a black hole of ~ 20 ± 5 Mq was formed. 

The above evolutionary history is only only one of the 
possible scenarios. As noted by Gies (2004, private commu- 
nication), the above scenario implies that about half of the 
total initial mass of the binary system was lost in the form 
of stellar wind. This implies, in turn, that the initial or- 
bital period was < 3 days, which means that the system 
was very tight. One cannot exclude that the initial masses 
were smaller, the initial orbital period longer and the system 
passed through the mass transfer and the common envelope 
phase. 



7 CONCLUSIONS 

(1) The calculations modelling the evolution of HDE 
226868, under different assumptions about the stellar 
wind mass loss rate, provide robust limits on the present 
mass of the star. For the most likely intervals of the 
values of the distance and of the effective temperature: 
1.95 to 2.35 kpc and 30000 to 31000 K, the mass of 
HDE22868 is 40 ±5 Mq. Extending the intervals of 
these parameters to 1.8 to 2.35 kpc and 28000 to 32000 
K, one obtains the mass of the star in the range 29 -r- 50 

Mq. 

(2) Including the additional constraints resulting from 
the observed properties of the binary system HDE 
226868/Cyg X-1, one can estimate the mass of the black 
hole component. For the most likely values of the pa- 
rameters mentioned in the item (1) above this mass is 
20 ±5 Mq . For the extended intervals of the parameters 
the mass is in the range of ~ 13.5 ~ 29 Mq. 

(3) The distance to the binary system has to be in the 
ranges 1.80-M.88, 1.98^2.15, 2.04-^2.22, 2.074-2.27 and 
2.21-^2.35 kpc, for the effective temperature of HDE 
226868 equal to 28, 30, 30.7, 31 and 32 x 10^ K, corre- 
spondingly. 
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Table 2. Parameters of the selected evolutionary models of the binary system HDE 226868/Cyg X-1. 
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Te Mopt d Mo fsw Ropt logL -M /rl * 

[lO^K] [Mo] [kpc] [Ms] [Rs] [Lo] [10-« [°] [Mq] 

Al . /yr] 



Observational paraiiiotcrs 



22.77 5.617 2.60 > 0.95 33 

±2.3 ±0.13 -1.3,+2.6 ±5 



Acceptable models 



28 
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1.80 


33.5 
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5.330 


1.84 


0.95 


34.6 


13.6 
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19.9 


28 


30.31 


1.88 


35 
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20.53 


5.367 


2.12 


0.98 


33.9 


14.2 
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(1) Mo denotes the initial (ZAMS) mass of the optical component, /gw denotes the multiplying factor 
applied to HPT formula; other symbols have their usual meanings 
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